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■ We address the problem of potential non- universality of the leading 1/Q power 

corrections to jet shapes emerging from the non-inclusive character of these ob- 
. servables. We consider the thrust distribution as an example and analyse the 

non-inclusive contributions which emerge at the two-loop level. Although for- 
mally subleading in a s , they modify the existing naive one- loop result for the 
expected magnitude of the power term by a factor of order unity. Such a pro- 
motion of a subleading correction into a numerical factor is natural since the 
non-perturbative power terms are explicitly proportional to powers of the QCD 
>• ! scale A which can be fixed precisely only at the two-loop level. The "jet-shape 

scaling factor" depends on the observable but remains perturbatively calculable. 
Therefore it does not undermine the universal nature of 1/Q power corrections, 
which remain expressible in terms of the universal running coupling and universal 
soft-gluon emission. 
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1 Introduction 



The analysis of power-behaved contributions to infrared and collinear safe characteristics of 
hard QCD processes has recently been developed as a method to quantify non-perturbative 
effects It has been recognised that the leading power corrections to infrared/collinear 

safe quantities are determined by an incomplete cancellation of soft-gluon contributions. 
A QCD process at a hard scale Q ^> A is determined by the small space-time intervals 
between the quarks and gluons involved. In spite of that, at the level of power suppressed 
corrections, the cross section of such a process may acquire a contribution from a soft gluon 
which travels a finite distance lfm ^> 1/Q and thus is sensitive to the non-perturbative 
interaction domain. 

The sensitivity of a given observable to soft-gluon radiation can be studied perturbatively 
to determine the power n of the expected leading Q~ n correction. Such theoretical input 
is valuable on its own, especially for jet physics which deals with intrinsically Minkowskian 
quantities and therefore typically receives no hints from the Euclidean-based OPE technology 
0. For example, an analysis of large-distance soft-gluon effects predicts that the properties 
of jets assembled with the use of the (once standard) JADE jet finder are contaminated 
by confinement effects at the level of Q^ 1 corrections, while the Durham (k^_) algorithm 
produces jets whose rates will depart from perturbative predictions only at the Q~ 2 level [Q. 

More challenging is the prediction of the magnitude of expected power corrections. For a 
single observable this would be impossible given the present state of the art. Nevertheless it is 
natural to expect that soft-gluon universality implies proportionality between the magnitudes 
of power-behaved terms in different observables ||. Along these lines one aims to predict the 
relative sizes of power corrections to observables belonging to the same class, linking together, 
for example, jet mass and thrust distributions with 1/Mq effects in heavy quark energy 
spectra MM (n = 1) or power corrections to the Drell-Yan .fT-factor |||| with those in the 
Gross-Llewellyn-Smith DIS sum rule and e + e~ fragmentation functions (n — 2). 

However, some doubt has been expressed as to the universality of power corrections to 
jet shapes because the latter are not fully inclusive with respect to final-state branching of 
the soft gluon under focus. To the best of our knowledge this question was first considered 
by Nason and Seymour in ||. 

The notion of soft-gluon universality contains two ingredients. First, it exploits the 
universal character of the soft-radiation matrix element in the spirit of Low's theorem [[10 . 
Secondly, it requires universality of the gluon radiation intensity (running coupling). 

The running QCD coupling in hard processes emerges after inclusive integration over the 
gluon decay products, resulting typically in ac B (k±) with k± the gluon transverse momen- 
tum [PH . In QED it suffices to consider the virtual photon self-energy blobs ("bubbles"), 
while in the QCD context there are also the corrections to the gluon emission vertex. The 
dispersive representation has been used in Jf] to give meaning to the QCD running coupling 



at low scales PJl^j. 

The value of a jet-shape variable depends, however, on the kinematics of all final state 
particles (hadrons, partons). As long as the observable is collinear safe the (quasi-)collinear 
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gluon decays, which are the main contributors to the running a s , do not affect the observable 
and therefore can be treated inclusively. However, in the "large-angle" gluon decay region, 
the inclusive interplay between real and virtual gluon loops gets broken because a real decay 
may correspond to a value of the observable different from that for the case of a virtual loop. 
These contributions seem to undermine the above-mentioned universality. 

In the case of the thrust distribution, for example, the naive inclusive treatment applies 
as long as the decay products fall into the same hemisphere. Then the contribution to thrust, 
T, can be attributed to the parent gluon, and integration over the virtuality of the latter 
results in the running coupling in a standard way. 

The existence of the part of the decay phase space in which the two gluon offspring 
partons (quarks or gluons) happen to fly into opposite hemispheres, has two consequences. 
First, this region has to be treated separately as its contribution to T depends on details 
of the gluon decay matrix element. At the same time the effective intensity of emission of 
the parent gluon as a whole gets modified. This is due to the fact that to reconstruct the 
running coupling one sums inclusively over the gluon decay products, while in this case one 
modifies the final state phase space. Thus the running coupling (its "spectral density" ) loses 
part of its standard support. 

It has been argued that such a non-inclusive contribution is higher order in a 8 . Its 
quark- ant iquark ("Abelian") part was studied numerically in || and found to be small. 

In this paper we analyse the question of (non-)universality of the leading power corrections 
to jet shapes and thrust in particular. We find that the naive inclusive treatment has to be 
modified by two effects. The first is due to non-inclusive configurations of the parent gluon 
decay (non- inclusive correction). The second is due to an incomplete compensation between 
real and virtual logarithmic terms coming from the soft-gluon splitting into gluons (inclusive 
correction). This correction is absent in the contribution from the gluon splitting into a 
quark- ant iquark pair. We show that both non-inclusive and inclusive corrections, although 
formally sub leading in a s , get promoted to a factor of order unity in the magnitude of the 
1/Q power correction. As a consequence, previously obtained results based on the naive 
inclusive treatment of an observable j acquire a "rescaling factor" (l + r^) which is calculable 
but depends on the observable under consideration. 

In this sense one can still consider power terms in jet observables as universal, i.e. based 
on the universal coupling and the universal soft-gluon radiation pattern, although their 
relative magnitudes differ from the naive expectations based on one-loop off-shell soft-gluon 
matrix elements. 

For thrust values close to 1, only soft-gluon radiation is essential. For T well below 1, one 
has contributions both from multiple soft-gluon emission and from rare iV-jet configurations 
(N > 2). The first contribution has been found || to be dominant at least down to T = 2/3, 
the minimum T for any three-jet configuration. As far as power corrections are concerned, it 
was shown in [§] that the order a s (Q 2 ) three-jet configuration does not develop a 1/Q power 
correction. 



Therefore to analyse the thrust distribution at two loops we take into consideration [ 1 3 



particle ensembles consisting of the primary quark- ant iquark pair accompanied by partons 
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originating from multiple emission of primary soft gluons and employ the existing all-order 
resummed perturbative prediction |L4 . 



We extract the leading 1/Q power term in the radiator, namely the exponent of the 
Mellin-transformed T-distribution. We do this in the framework of the treatment of non- 
perturbative power corrections to infrared/collinear safe observables introduced in Q], based 
on the notion of an infrared-finite strong coupling which differs from the perturbative form 
in the infrared region. 

At the two-loop level, non-inclusiveness of the thrust reveals itself and the characteristic 
QCD scale A becomes precisely defined through the fixing of the scheme for the running 



coupling. To this end we choose the CMW renormalisation scheme [I3|, in which a s is 
defined as the intensity of soft-gluon radiation. 

Section 2 contains the necessary ingredients for the two-loop analysis of the radiator of 
the thrust distribution: kinematics, structure of the soft-gluon emission followed by gluon 
splitting into qq and gluon pairs, scheme fixing, and thrust resummation. In section 3 
we recall the dispersive method of and its application to the thrust distribution in the 
naive inclusive approach JJ. Section 4 is devoted to the exact treatment of the problem, 
the analysis and calculation of the inclusive and non-inclusive contributions to the thrust 
rescaling factor (1 + r T ). Section 5 contains a discussion of the results and the prospects for 
the future. 



2 Two-loop analysis of thrust distribution 



2.1 Thrust: kinematics 



Consider a QCD hard process in which the final state consists of a leading quark- ant iquark 
pair with momenta p and p together with a system of N partons (quarks and gluons) with 
momenta k\, ■ ■ • , fcjv- It is convenient to cast the cross section in terms of the Sudakov (light- 
cone) variables defined with respect to the thrust axis of the event. To this end we introduce 
two light-like vectors 

pp ) p> . p2 = p2 = Q . pn + p> = g M ^ 2 (PP) = Q 2 = 1 , 
and represent the parton 4-momenta as 

K = fcP* + + (k l± r, aA = k 2 ± 

The leading quark momenta become 

p = (i-EA- Pp)p + + p± , 

1=1 

p = (1 - a i - a v)P + PpP + P± ; 



i=l 
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Adding the normalised longitudinal momenta of the secondary partons, k^j Q — | ctj — /3{ | /2, 
with those of the two quarks, for the thrust value we obtain 

T = ^2\ \ ai - f3i\ + -^2 (3i - (3 P ) - a p ] + -^ai - a p ) - f3 p ] 

i i=i i=i ^ 

= 1 - min{aj, A} - (a p + (3 P ) . 

i 

This expression implies that p and p momenta belong to opposite hemispheres, which is the 
dominant configuration (the correction is numerically small and relatively suppressed at least 
as (1 — T) when T is large). 

The values of the "small components" of the quark momenta, a p and f3 p depend on quark 
transverse momenta with respect to the thrust axis through the on-mass-shell conditions, 
p 2 = p 2 = 0. Transverse quark recoils are typically of order 

P 2 ± ~ Pi ~ Vn ■ k\ , 

with k±_ characteristic transverse momentum of a secondary parton. We have kept in the 
estimate a (logarithmic) enhancement factor due to the parton multiplicity N. 

As we shall see the leading power correction comes from the region ctj ~ of the order 
of the transverse momentum (in units of Q). Therefore in ( |2.1| ) we shall neglect a p and /3p 
which are of the order of the transverse momentum squared. 



2.2 Gluon splitting 

Two-body variables and phase space. The two-parton phase space is 

Introducing the relative fraction of a-components, z, and the invariant mass squared of the 
two-parton system m 2 we write: 

a = a% + a 2 , Oil — zoc , a 2 = (1 — z)a , 
m 2 = (kx + k 2 ) 2 = (ai/3 2 + a 2 (3i) - 2k 1± ■ k 2± = z(l - z) 



k 



1± K 2 ± 



Z 1 — Z 



In terms of the scaled parton transverse momenta, 



ku_ k 2± _ _ 

<h = — > Q2 = - — , g = qi-q2, 

z l — z 

the phase space takes the form 

dT 2 (h, k 2 ) = ^dz^^ dm 2 6 ( - & - q 2 f) . (2.2) 
a IT 71 \z{l—z) J 

For later use we also define 

k± = k 1± + k 2± = zqi + (1 - z)q 2 ; 



k 2 + m 2 „ z 2 l — z 2 (2-3) 



£ = + & = — , A = -q( ; A = — q z 2 ■ 

a a a 
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Matrix elements. The squared matrix elements for the production of two gluons or a 
quark pair with 4-momenta k±, k 2 read, in the soft limit {Qkj) <C Q 2 (see ||13|| ) 



M 2 (k u k 2 ) = 4C 2 W X W 2 + 4CV (M 2 g + M t 
A(M 2 gg + M 2 qq ) = C A (2S + H g ) + n f H q . 



(2.4) 



The first term, proportional to Cp, describes the independent (Abelian) emission of two soft 
gluons off the hard qq, where 



pp 2 
(pk i )(k i p) kj ± 



is the standard dipole radiation cross section. 



The second (non- Abelian) term, proportional to Cp Ca, contains two pieces. The S term 
describes the leading infrared contributions which is singular in the ratio of gluon energies: 

S EE W W + W 2 1 ~WxW 2 , 

where 

w _ (PP) = 1 

13 (pki)(kikj)(kjp) m 2 a^j 

stands for the ordered 2-dipole emission. 

H g and H q are responsible for parton configurations with energies of the same order. The 
"hard" parts of the two-gluon (H g ) and the two-quark (H q ) contributions read 

H g ee SR + J 2 -A^\, (2.5a) 

H ° - - }2+2 v§y (2 ' 5b) 



Here 



J 2 



R 



\ 2 / 

(ph)(k 2 p) - {pk 2 )(kip)\ ( ai(3 2 - a 2 (3i 



\ (k\k 2 )(pk)(kp) J \ m 2 a[3 
ipkx){k 2 p) + (pk 2 )(kip) ax/3 2 + a 2 Pi 



(2.6) 



(pk)(kp) a(3 
and W{k) = (pp)/(pk)(kp) is the dipole emission of k = k\ + k 2 . 



We may write the squared matrix elements in the following form: 

M 2 +M 2 = -, = M rsp-, (2.7) 

99 qq m 2 (zql + (1 -z)ql) V ; 

where 



M 2 = 2C A (2S + H g ) +2n f H q . (2.8) 
The distributions determining M 2 are functions of the dimensionless ratios u 2 with 



u = — 

\q\ 



Namely, 



2S 

Hg 

H a 



zu\ + (1 — z)u\ ( 1 
z(l-z) 

-2 + z(l-z) 



z(l-z) 



K - ul? 

Zuj + (1 — z)u\ 
ZU\ + (1 — Z)u\ ' 



2 2 1' 

ufuij 

u\ + u\ 



The collinear limit (q 2 oc m? — > 0) corresponds to taking w 2 



(2.9a) 
(2.9b) 
(2.9c) 

oo while keeping fixed 



2 ' 1 2 2 

Ml «2 U{U2 



Uo 



the modulus of the vector difference, 



(ttl - M 2 ) Z 



see ( p.2|) . Observing (see Appendix A) that the azimuthal average of the ratio 

-> 2, 



(«i - 



we derive the parton splitting functions 

2C A (2S + H g ) 
2n f H q 



AC, 



2n/(l -2*(1 -z)) . 



2 + z(l 



(2.10a) 
(2.10b) 



Note that the latter expression is twice the standard g — > qq splitting probability because 
thrust as an inclusive observable does not differentiate between quarks and antiquarks, and 
the symmetry factor 1/2! will be included in the thrust trigger-function, see ( |2.23|) below. 



2.3 Gluon decay and the /^-function 

The /3-function is related to the integrated two-parton production probability. Keeping fixed 
the invariant mass and total transverse momentum of the pair, and integrating M 2 over 
the azimuthal angle with respect to k± and over the relative momentum fraction z, it is 
straightforward to derive (see Appendix A) 

f 1 dz (** — (M 2 + ML) = „„ 9 1 ^ 

2! Jo Jo 2vr V *g m) ^(kl + m 2 ) 

where 

The logarithmic singularity at m 2 — > comes from the soft singularity in the g — > gg matrix 
element (splitting function) and is cancelled by the corresponding virtual correction to the 
one-gluon radiation vertex. The latter can be written as an integral over an intermediate 
gluon virtuality /i 2 , 



u 0+2CAl jm±^i\ , (2 . n) 
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Here we have chosen the form of the logarithmic term so as to cancel exactly the corre- 
sponding real emission contribution in ( |2.11| ). The mismatch function s(x) vanishes at the 
origin, s(0) = 0, and depends on the scheme one adopts for calculating the divergent vertex 
correction. 



2.4 Fixing the scheme for a s 

To fix the scheme for defining the perturbative expansion parameter a s , we have to consider 
first the relation between particle production at order a 2 an d the coupling. To this end we 
address the non-Abelian contribution to an inclusive quantity such as the quark anomalous 
dimension. 

The inclusive probability for producing a soft gluon or a two-parton system with gluon 
quantum numbers reads 

da dk 2 
a k 2 

where the "anomalous dimension" 7 becomes 



dw = 4C F — ^7(0, (2.14) 



' - 2+**+r {-^^m^i . ( , 15) 



The singular logarithmic term cancels against that in x{k±), and we arrive at 

Here c s is a constant given by the integral of the scheme-dependent piece of Hu, see (|2.13| ): 

f°° dx 
Cs = / -7—. — r s[x) . 
Jo x(l+x) 

First we observe that, to the necessary accuracy, 7 satisfies, 



2 d ( cXsikp V _ 2 d a s (k 2 ±) 



This means that we can choose 7 to represent the running coupling a s . Equation ( p.16 ) then 



becomes a dispersive relation for a s (k"j_), in which the combination of «o and c s a 2 defines 
the coupling constant at k\ = 0: 

ajk 2 :) If., r°° dm 2 k 2 , . 9 . 1 1 f°° dm 2 , 9n 



47r 47r [ s jo m 2 (/c^_ + m 2 ) s j 47r 7o /c 2 + m 2 

(2.17) 

where p s ijn 2 ) is the spectral density, to be discussed later. 



This corresponds to choosing the CMW renormalisation scheme JL5| in which a s is defined 
as the intensity of soft-gluon radiation. In another scheme, e.g. MS, the anomalous dimension 
in the soft limit would contain both a s and a 2 terms. 
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2.5 Thrust: resummation 



For a qq + iV-parton final state the thrust distribution becomes 




N 

dT N M 2 ( Ql , q 2 ,... q N ) 5(1 — T — mm{a u ft}) , (2.18) 



with M the iV-parton matrix element and dT^ the corresponding iV-body phase space. In 
the soft approximation (1 — T 1) the final partons are either relatively soft gluons emitted 
off the quark-antiquark line (primary gluons) or their decay products. The matrix element 
for an ensemble of primary gluons is very simple in the soft approximation, as it is given 
by the product of independent radiation amplitudes (essentially an Abelian pattern). By 
taking a Mellin representation for the ^-function, 

( N \ t dv N 

8 I 1 - T - minfaii, ft} J = J — exp{v(l - T)} J[ exp{-z/min{a;, ft}} , 

the thrust distribution takes the form 

At next-to-leading order, the "radiator" contains two terms, TZ(v) — TZ\ + 7Z 2 , due respec- 
tively to one- and two-parton contributions. 



One-gluon contribution. The one-gluon contribution, TZi, is given by the following for- 
mal expression: 

where the B-function selects the smaller longitudinal component (a < (3 = kj_/a). The term 
with exp(— vol) corresponds to the case in which the gluon is emitted, while —1 corresponds 
to the case in which the gluon is in a virtual loop. The function x{k±) * s the order a 2 con- 
tribution that describes the /cj_-dependent one-loop vertex correction to one-gluon emission, 
given in ( |2.13| ). 

Using the expression for the virtual correction, ( |2.13| ), we can split IZi into two contri- 
butions 

K x = Tin + Hvi , 

where 

Kll = SC F /' * I d -$ - l) 9(k ± -a)-^f-, (2.21a) 

The weird factor a s (0) in (|2.21a|) , which stands for the two- loop value of the on-mass-shell 
coupling, is there simply because a real gluon has been emitted. The contribution 7Zi 2 is 
also ill-defined. 



8 



As we have seen above, IZn, is cancelled by the contribution from two-parton gluon 
splitting in the collinear limit in (|2.11 ), resulting in a finite answer in which a s (0) gets 
replaced by the running a s . 

The second contribution, TZu, taken together with the logarithmic piece TZ 2 2 of the two- 
parton emission in TZ 2 , defined below, produces a finite inclusive correction to the answer. 



Two-parton emission contribution. The second part of the radiator describes gluon 
decay into two quarks or into two gluons. The latter is described by the non-Abelian part 
of the qqgg matrix element M 2 g of (2A), proportional to CfCa, see [13|]. We have 

n 2 = 4cv / dv 2 ) (Mg g (ki, h) + M 2 gq (h, h)) n T . (2.22) 

Here essential higher-order virtual corrections have been taken into account, leading to the 
running of the coupling with the invariant mass of the virtual gluon, \a s (—m 2 )\ 2 = a 2 (m 2 )(l+ 
Oin:)). 

The "thrust trigger-function" Qt in (|2.22|) determines the contribution to the thrust of 
any two-parton configuration and is given by 

= 2j [exp(-z/min{o;i,/3i} - z/ min{a 2 , ^2}) - 1 ] • (2.23) 

We note that having included the symmetry factor 1/2! for both gluons and quarks we have 
implicitly defined M 2 Jki, k 2 ) to describe the probability of finding either a quark or an 
antiquark with momentum k\. 



Splitting the answer. Now we are ready to assemble the pieces of the radiator into three 
contributions each of which is explicitly finite. 



The "thrust trigger-function" Qt in ( |2.22|) depends on the azimuthal angle between the 



partons, 0, and does not allow the integration in (|2.11| ) to be performed. It can however be 
done, in a naive ("inclusive") approximation in which one employs a simplified version of 
the IMactor, namely, 

4 simp) = 1 [exp(-z/min{a, (3}) - 1 ] , (2.24) 

where a = a.\ + ct 2 and (3 = fix + /3 2 . This corresponds to the contribution to the thrust 
from the parent gluon, disregarding the possibility that its offspring may go into opposite 
hemispheres. The regions a > (3 and [3 > a contribute equally, so that it is sufficient to 
consider only the latter: 



Using q and k± as integration variables, we obtain 
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with the angle between q and k±. 



The integral in ( |2.26| ) has been calculated above in (|2.11|) . It consists of two terms 



n { r mp) = n 21 + n 22 : 



72 



^ ^/"(^-O/^g^) ■(-« m 2 + (2-27.) 

Thus we may represent the answer as a sum of three finite terms: 

n = k {0) + n {i) + n in) , (2.28) 

where 

= ftn + 722i , (2.29a) 

72 (i) = 72i 2 + 72 22 , (2.29b) 

km = n 2 ~n { 2 simp) . (2.29c) 



In the first contribution ( |2.29a| ) 72n involves a s (0) which is cancelled by 7Z 2 \ to give the 



standard "naive" answer, as we shall see in the next section. The second contribution, 
we shall refer to as an "inclusive correction". Here both 72i2 and 1Z 22 contain m 2 — > 
logarithmic singularities which, however, are damped by the difference of the two simplified 
"thrust trigger functions" 

4 simp) - 4 simp) (m 2 = 0) oc Q(k 2 ± + m 2 - a 2 ) - Q(k\ - a 2 ) . 

Similarly, the non-inclusive correction (|2.29c|) is also finite as it involves the difference of the 
true and the simplified trigger functions, 

5n T = fir - fij imp) , (2.30) 

which also vanishes as m 2 —> 0. 



Now we are ready to recall the dispersive method and apply it to the calculation of ( 2.29|) , 
starting from the main term 72 (0) . 



3 Naive power correction for T in the dispersive method 
3.1 Dispersive method (recollection) 

In the approach of @ the dispersive representation for a s is implemented, in which the 
dispersive variable m 2 acts as a "gluon mass" in the matrix elements (see |2.17| ): 

a ' (k2) - r *" 2 rA K)=r^4,Mm 2 ). (3.1) 



k 2 Jo k 2 (k 2 + m 2 ) Jo (m 2 + k 2 ) 2 
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At one-loop level one has 



ps(m 2 ) = d a eS (m 2 ) = -^a 2 s {m 2 ) + ... (3.2) 
a In m z An 



The perturbative answer, Oj(Q 2 ), for an observable j is then formulated in terms of the m 2 
integral 



with the "characteristic function" one derives by computing the m 2 -dependent soft-gluon 
radiation matrix element. For an infrared/collinear safe observable the m 2 integral converges 
and is mainly determined by the region m 2 ~ Q 2 , the typical hard scale of the process. This 
region reproduces the standard one-loop perturbative prediction 

Of T = a eS (Q 2 )F(l) « a s {Q 2 )F{l). 

The genuine non-perturbative component of the answer, 50^^ , is triggered by the non- 
perturbative component of the effective coupling 5a e ff(m 2 ). Imposing the ITEP/OPE re- 
striction H that a s (k 2 ) in the ultraviolet region receive no power correction (k 2 )~ p due to 
the large- distance interaction domain, one derives from ( |3.1| ) that the integer moments of 
the non-perturbative effective coupling vanish: 

roc dfji^ 

— -Sa eS (m 2 ) (m 2 Y = , i = 1, 2, . . .p . 
m z 

Thus, the power corrections are obtained from the non-analytic m 2 — > behaviour of the 
one-loop matrix element (characteristic function). For example, the behaviour 50^^ ~ Q^ 1 
in the thrust distribution originates from 

T T [m 2 IQ 2 ) - 2 -^M, m 2 -0. (3.4) 

7T V Q 

The magnitude of the power correction is then expressed in terms of the m-moments of an 
"effective charge" 5a e g(m 2 ) describing the intensity of non-perturbative gluon radiation: 

SOf F (Q 2 ) = r ^5a eS (m 2 ) t 3 {m 2 IQ 2 ) z> , (3.5) 

where 

A 2p . q = ^ T ^ 5a eS (m 2 ) • (m 2 f In* m 2 (3.6) 
2tc Jo m z 

(for p half-integer or q > 0). In particular, for the case of thrust from ( |3.4j ) we deduce: p — 5, 
g = 0, and the observable-dependent constant Cj is = 4. 



3.2 Inclusive treatment of the thrust distribution 



In an inclusive approximation one employs the simplified version of the Q- factor ( 2.24 ). The 
contribution 7£ 2 i ( J2.27a[ ), proportional to /3 , is the only one which has been dealt with in 
the past [ 



Making use of (|3.2|) it can be written as 



K 



2Cf r da 



21 



7T 



a 



(e~ ua - l) 



x -da eS (m 2 ) 



dm 



dm 2 



dk 



2 Q{k\ 



+ m 2 — a 2 ) 



k\ + m 2 



(3.7) 
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To reconstruct the known result we integrate by parts. Since a e ff(0) = a s (0), the surface 
term cancels TZn, given in ( |2.21a| ), and we obtain 

tt<°> =n u +n 21 = 2 -^ (e— - l) r a eff (m 2 ) 6(m 2 - a 2 ) . (3.8) 

The normal perturbative prediction is obtained by integrating over a between \ jv <C 1 and 
1. In order to examine the non-perturbative contribution, one substitutes for ct e s its non- 
perturbative part, 5a e f[, which is concentrated at small m? . This implies that a is small, 
allowing us to expand the exponential term, 

SR® = -J 1 ^- r^-m 5a eS (m 2 ) ee -v 4A lt0 , (3.9) 
7i Jo m 2 

resulting in a non-perturbative shift of the thrust distribution, proportional to the first 
moment of 5a e fr 0, 

da(T) _ daWjT-AT) (0) _ 4 A lfi 

dT ~ dT ' 1X1 ~ Q ■ {6A()) 

Taking into account subleading a 2 -effects, the naive inclusive prediction ( |3.10|) gets modified 
by the thrust rescaling factor tt' 

AT = i^(l + r T ). (3.11) 

This factor originates from the real and virtual two-loop matrix element and depends on the 
observable under consideration. 

In what follows we analyse and compute the thrust rescaling factor. In our framework, 
it consists of two components: an inclusive r^ and a non-inclusive r^ corrections 



4 Exact treatment 



4.1 Inclusive correction 



As we have seen above, the logarithmic terms in two-particle production 7Z 2 2 ( [2.27b|) and 



the vertex correction 72-12 ( ^-21 b| ) cancel in the fully inclusive anomalous dimension. 

In a less inclusive quantity, like thrust, the cancellation is no longer complete: diver- 
gences cancel but a finite correction remains. Indeed, in the real emission contribution the 
kinematical restriction involves the mass of the pair and reads 



6 yy k\ + m? — <\ 

In the virtual correction, however, we deal with a real gluon restricted simply by 

6 (\fkZ-a 
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As a result the finite correction which emerges reads 

r<« =i6c F c, I * - 1) / -I^L f ^) 

j a v ' 7 m 2 [kj_ + m 2 ) y 4/T y 

fcKfci + m 2 ) 



(4.1) 



e -e(k ± -a) 



Expanding the exponent and performing the a-integration we arrive at 

dm 2 ( a s (m 2 )\ 2 r°° dk\ 



J m z \ 47T y jo 



o i + m 2 



A; 2 +m 2 — k±) In 



fc 2 (fei + m 2 ) 



where we have extended the /c^-integration to infinity since it converges. Now we use 
and integrate by parts to arrive at 



C 'a 

8Cp— —v I dm 

Po 



,a eS (m 2 ) d 



An dm 2 



with 



,(0 



00 xdx ln[x 2 (l + x 2 



o 1 + x 2 x + \/l + x 2 
Substituting Sa e s for a c s we get the power correction 

(i) 2C F C A f dm 2 



{2mc (!) } 



3.299 . 



sn {t) 



-V ■ c 



m <5a e ff(m 2 ) . 



(4.2) 
(4.3) 

(4.4) 



Recalling the definition of the non-perturbative moment, ( |3.6| ), the relative correction to the 
shift in the thrust distribution ( |3.10|) becomes 



AT® 



4 ^i.o r « 







T J 



c 

r w = _d c (0 = L10 o (0.900) for n/ = 3 (0) . (4.5) 

Po 



Notice that the inclusive correction is of the same magnitude as the naive term. 



4.2 Non- inclusive correction 

Since the leading power correction originates from the region of large gluon radiation angles, 
a ~ (3 ~ m, k± < m, the kinematic region of offspring partons moving into opposite 
hemispheres will give an essential correction. 

This correction is taken into account by in ( |2.29c| ), given by the expression Q2.22| ) 
with the thrust trigger-function Qt replaced by 5Qt fl2.30|) . In order to obtain the first power 
correction it suffices to expand the exponents in ( |2.23| ) to first order as above. Straightforward 
calculation results in 

roo (^q/ _ / 

/ — 5Q ~ -vQt, Q T (q 1 ,q 2 ) = zq x + (1 - z)q 2 - J zq\ + (1 - z)q\ . (4.6) 
Jo a v 

The expression for IZ^ becomes 

7J(»> = -„ iC F J dm' 2 / ik l £ dz J d A (Ml + l£) S1 T . (4.7) 
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Now we use Q3.2|) and integrate by parts. Taking into account that at m = 

^t(<?i=<?2) = , 

the surface term vanishes and we get 



The expression in the curly brackets has dimensions of mass. Q 2 enters only in the upper 
limit of the k\ integration which is convergent in the ultra-violet region. Therefore we 
integrate up to infinity so that the result does not depend on Q and is given by 



{(PD) = 2r| n) m + 0{m 2 /Q) 



with a number still to be calculated. In the m — > limit this leads to 



5n (n) = _ v 4^0 r W (4 g) 



As a result, the non-perturbative shift in the thrust distribution (|3.11 ) acquires a contribution 



AT W = (41Q) 

In what follows we construct the relevant matrix element and compute the non-inclusive 
component of the thrust rescaling factor . 

4.3 Evaluation of r$ 

We are now in a position to discuss the determination of the non-inclusive thrust rescaling 
factor tt- From (|4.8|) , we have that 

( n)_i rd 2 k ± d0 M 2 nr(gi,g 2 ) uu , 

T ~mJ vr 2k zq 2 + (l-z)q 2 A(3 ' 1 j 

To express the integration measure in terms of u% and u 2 we write 

d 2 k ± d(f) d 2 k ± d 2 q 2 _ 2 d 2 q 1 d 2 q 2 2 ^ 2 
— = 8{q - (gi - g 2 ) ) = 5{q - {qi - q 2 ) ) 

71 Z7T 71 71 TX TX 

2 2 2 2 mdii! u 2 du 2 
= qidqi <hdq 2 - dtp 6{q -q 1 -q 2 + 2<?ig 2 cos 0) = , . 

7T 7T 2u 1 M 2 |Sin0| 

Using the 5-function condition to express sin we arrive at 

d 2 k± d<j) Aq 2 u\du\ u 2 du 2 



71 2-K 71 yfj 

where the Jacobian factor J(ui,u 2 ) is 

J — ({ui + u 2 ) 2 — 1) (1 — (ui — u 2 ) 2 ) . 



(4.12) 
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Making use of the fact that £1 in Q4.6|) satisfies 

~ ~ 771 ~ 

^r{q\,qi) = q^r(ux,u 2 ) = , Q T (ui,u 2 ) , (4.13) 

y/z(l-z) 

we finally obtain 

If 1 dz f°° f°° , 1 

:— - / / U X dU X / W 2 <™2 — Fr= 

^A) -/o . / z n - z ) Jo Jo vj 



T 



X ^ + (l-^| fir(Ul ' M2) 



(4.14) 



where the integrals run over the region J > 0. This region looks simpler in terms of u± 
Ui ± M 2 : 

dui du 2 i f 1 c?tt_ r 00 d« + 



V^7 ^" 2 7-1 /l _ M 2 h Jy2 _ 1 



(4.15) 



Convergence of r^p . It is necessary to check that the u + — > oo and z(l — 2;) — > regions 
do not lead to divergences of the integral. This is shown in Appendix B. 

Numerical integration results in the value 

4 n) = — (-1.227 CU + 0.365 CU - 0.052 n/) , (4.16) 

Po 

where the three terms originate respectively from the soft gluon, hard gluon and quark 
matrix elements (|2.9| ). This gives 

rP = -0.710 (-0.470) , for n f = 3 (0) . (4.17) 

We have that the non-inclusive correction alone would lower the "naive" expectation for the 
factor in the leading 1/Q power correction to thrust by about 30%. However, assembling 
the inclusive and non-inclusive corrections to the thrust rescaling factor we finally obtain 

r T = r f +4 n) = i- [1.575 C A - 0.104 n/] = 0.490 (0.430) for n f = 3 (0) . (4.18) 

Po 

The conclusion is that the thrust rescaling factor 1 + tt increases the naive value by 50%. 
The most important correction is the inclusive one ( |4.5| ). 



5 Discussion and conclusions 



Within the dispersive approach of power corrections to perturbative QCD predictions, 
and to jet-shape observables in particular, are expressed in terms of moments of the non- 
perturbative component of the effective coupling, 5a e ^{m 2 ). 

To fix the absolute normalisation for the leading power corrections to jet-shape observ- 
ables a two- loop analysis has to be carried out: since power corrections are proportional to 
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powers of the QCD scale A, it is necessary to fix the latter to control the magnitude of the 
non-perturbative correction. 



In this paper we have examined the thrust distribution in the high-thrust region to two- 
loop accuracy to extract the 1/Q power contribution. To do this it suffices to consider 
radiation of soft gluons followed by their two-parton decays. 

To this end we have used the physical CMW scheme |H| in which a s is defined as the 
strength of inclusive soft-gluon radiation (the magnitude of the singular part of the two-loop 
quark anomalous dimension). 

Given AT^ as the result of the naive inclusive one-loop treatment, we have found that 
the rescaling factor 1 + ry for the thrust power correction, AT = (1 + r^)AT(°^ is rather 
large. This large correction comes mainly from the gluon splitting into two gluons. 

The correction tt originates from two sources. A positive term (the "inclusive correc- 
tion") comes from an incomplete compensation of the logarithmic contributions in soft-gluon 
splitting, due to the kinematical difference between virtual (massless) and real (massive) 
gluon radiation ( |4.1| ). This gives the most important correction (see |4.5| ). A negative term 
rj^ (the "non- inclusive correction") is due to the kinematical region of the 4-parton phase 
space in which the gluon decay products fly into opposite hemispheres (^4.7| ). The contribu- 
tion to this last term coming from the decay of the primary gluon into a quark-antiquark 
pair proves to be positive and numerically very small. This agrees with the result of the 
analysis of Nason and Seymour || . 

It is the region of large gluon emission angle, 9 ~ 7r/2, that triggers the non-analyticity 
in the gluon virtuality m 2 and thus power corrections. This makes opposite-hemisphere con- 
figurations of the decay products quite common. But it so happens that the most important 
correction is the inclusive one. 

Since the moments of 5a e s(m 2 ) are not calculable given the present state of the art, 
they must be determined phenomenologically. Therefore to be of any practical value, the 
programme of calculating power corrections at two loops must be extended to other ob- 
servables driven by the same (first) moment as the thrust. In a forthcoming publication 
we will generalise the present analysis to energy-energy correlations and jet-broadening in 
e + e~-annihilation JTB| . 



For a general observable j we will have a rescaling factor 1 + rj, where rj originates 
from the primary soft-gluon emission and its two-parton decay convoluted with the trigger 
function specific to the observable j. Since rj depends on the observable, one could argue 
that universality of the leading power correction is violated. On the other hand, we would 
prefer to consider universality as still being valid. Indeed, the ingredients of universality, 
which are the universality of soft-gluon radiation and universality of the coupling, remain 
intact in the two-loop analysis. Moreover, only at the two-loop level can universality be 
given a precise meaning, since only at this level can one define the QCD scale and thus 
the coupling. The fact that the rescaling factor is observable dependent does not break 
universality, since rj remains under perturbative control. 

Let us stress that the two-loop analysis has led to a term in the rescaling factor which 
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does not contain any small parameter and is just a number. One may wonder whether 
higher- loop calculations will give corrections of the same order as r r We argue that this is 
not the case and that terms higher than two-loops give "subleading" contributions, i.e. with 
additional powers of the perturbative coupling. 

To see this observe that the (£+2)-loop matrix element would lead to the appearance of an 
m 2 -integral of af +2 (m 2 ). By using ( |3.2j) and performing the analysis as before, one finds that 
the corresponding contribution to the non-perturbative power correction will be proportional 
to al s 5a e ff and thus will be down, with respect to the two- loop result, by an extra power I 
of the coupling. Strictly speaking, such a perturbative argument is far from being perfect 
since the coupling here enters at a small scale. Whether the higher order contributions are 
actually small depends on the numerical value of the "perturbative" effective coupling a c s/n 
at low scales. The large value of the two-loop correction r T indicates that 5a e ^ is smaller 
than was indicated by the naive treatment. 
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of this paper, discovered subsequent to an independent calculation of the n^-dependent piece 
in pi. 



A Two-parton production cross section, ( |2.11| ) 

Using the reduced squared matrix element M. 2 in (|2.7Q , we can write 

(a.i; 

where the distributions 2S, TC g and 7i q are given in (|2.9|) as functions of Uj = qi/q with 



<fi = k±_ + (1 - z)q, q 2 = k± - zq. 
Performing the integration over 0, the angle between k± and q, we obtain 

^ dd> 1 / 1 \ 1 



o 2vrg 2 \qf / \k 2 ± -(l-z) 2 q 2 



2n — 1 1 / 1 - Z Z \ (A 2) 

o 2tt q\q\ m 2 + k 2 _\ q\ q\j 

^dct>q 2 lq 2 \ 2 

O 2 = 1 + ( — / B \ Z ) ' 

o 27cq£ \qi z 



where 



b(z) = e(,v - ki) = ©^-p^: 
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Similar expressions hold for (1/gf) and (q%/q 2 )- By using (|2.9|) we find 



1 (25) l-B(z)-B(l-z) 



2! v ' z(l-z) 



l //1/x rt „ N 5(2) 5(1 -z) m 2 l-6z(l-z) 
i<K,) = i-,(l-z) m> 1 - 8 »( 1 -*) 



2! 2 v y + m 2 2 

The last two terms in (Tig) and (Tig) vanish upon ^-integration and we find 

f X dz f 27T ^ {C A (2S + H g ) + n f H q } = 2 I' dz{c A {-2 + z{l~z)) + r ^{l-2z{l-z)) 
Jo Jo In Jo L I 



= -/5 + 2C A ln 
This gives ( |2.11 



\-B{\-z)- \B{zj) +^(l- - - *)) } 
fc 2 (fcj. + m 2 ) 



m 4 



(A.3) 



B Convergence of 



To check that the integral ( |4.14| ) determining converges, we first consider the large- 



it., 



behaviour of the integrand. The reduced matrix element ( j2.8|) has a finite limit at u + = 00. 
The integration measure in ( |4.15| ), together with the ratio UiU2/{zu\ + (1 — z)v%), behaves 
like du + /u + . The function Q that triggers the non-inclusive correction to thrust, ( |4.6|) , can 

be cast as 

Q(ut, u 2 ) = \ lu+ + (2z - l)u- - ^u 2 + + u 2 _ + 2u + u_(2z - 1)| 

2z(l-z)u 2 _ (B-l) 



(B-2) 



u+ + (2z - 1>_ + + (2z - l)u_) 2 + 4z(l - z)ul 

The crucial point is that it vanishes for large u+, 

~ u + ^ao z(l - Z)U 2 _ 

il{u 1 ,u 2 ) ~ , 

thus ensuring convergence. 

By inspecting (p.2|) we conclude that Cl also vanishes at z(l — z) = since 

nr ^ - z ) u - ,-a q \ 

fi(Ml ' M2) * , + + (2, - 1),, • (R3) 

In general, this is sufficient to compensate for the soft singularity of the matrix elements 

M 2 oc 1 



z{l-z) 
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Special care should be taken when analysing the edge of the phase space which corresponds 
to one of the partons being collinear with the quark direction: 



u + — 1 ~ 1 — \u~\ ~ J z{\ — z) <C 1 



where (|B.3|) vanishes only as y z(l — z). 

For z — > the potentially singular region is u± — > 1, which can be explored by writing 
u + = I + (p cos ip) 2 , U- = 1 — (psinip) 2 . 
For z and p small we have 

zu\ + (1 - z)u\ = ^p 4 + z + O (zp 2 ) . 

This combination appears in the denominators of the phase space and thus generates singu- 
larities for p — > and z — > 0. 

Consider first the region p 2 > y/z. Here the integrand of rjf 1 has the form 

dl ~ dz ^ 6 (p 2 - y/t] # M 2 . (B.4) 

For the soft reduced matrix element S one has 

2p 2 (cos 2 V - sin 2 V) + (p 4 ) (p 4 ) 



S 

z z 

where the first term is cancelled upon ip integration. The integration over p is convergent 
for p — > and gives a contribution with the integrable singularity 1/y/z. 

For the hard reduced matrix elements TC g and 7i q the dominant contribution in this region 
is a constant. Thus the p integration gives a contribution regular in z. 

Finally, the integration over p in the region p 2 < \fz leads to a contribution which, after 
angular integration, is regular in z. 
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